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ABSTRACT 

We present an approach for analysing the morphology and physical properties of Hi features near giant OB asso- 
ciations in M33, in the context of a model whereby the Hi excess arises from photodissociation of the molecular gas 
in remnants of the parent Giant Molecular Clouds (GMCs). Examples are presented here in the environs of NGC604 
and CPSDPZ204, two prominent Hn regions in M33. These are the first results of a detailed analysis of the environs 
of a large number of OB associations in that galaxy. We present evidence for "diffusion" of the far-UV radiation 
from the OB association through a clumpy remnant GMC, and show further that enhanced CO (1-0) emission appears 
preferentially associated with GMCs of higher volume density. 

Subject headings: galaxies: individual (M33) — galaxies: ISM — ISM: clouds — ISM: molecules — 
HII regions — ultraviolet: galaxies 



1. INTRODUCTION 

It is widely accepted that stars form from molecular 
gas in galaxies, and that the largest stars form in Giant 
Molecular Clouds (GMCs). Such massive stars produce 
prodigious quantities of far-UV (FUV) radiation, that 
will react back on the parent GMCs producing photodis- 
sociation regions (PDRs) on the cloud surfaces. Highly 
luminous PDRs very near to OB associations are a rich 
source of atomic, molecular, and solid grain spectral fea- 
tures, especially in the IR, and several prominent cases 
(e.g. the Orion region) in the Galaxy have been exten- 
sively studied both t heoretically and obseryationa lly (see 
e.g. the review by iHollenbach fc Tielensl Il999h . Less 
well appreciated is that Hi atoms are also produced in 
such PDRs in quantiti es sufficient to be d etectable in 21- 
cm radio emission (see lAllen et a l. (2004) and references 
cited there) , and that this can occur in situations of much 
more modest volume densities and far-UV fluxes than 
is usually required to produce detectable IR line emis- 
sion. Furthermore, such relatively low-density PDRs are 
physically large, rendering them more easily detectable 
by radio astronomy. 

In this paper, we present the first results of a study 
of Hi features in the the environs of OB associations in 
M33 (NGC 598) in terms of this PDR model. We ex- 
plain the observations in a consistent, simple physical 
context, u sing a method mo tivated by the initial dis- 
coveries of I Allen et al.l (|1986h . Those authors presented 
evidence that large-scale spiral arm features in the Hi dis- 
tribution of M83 (NGC 5236), offset from th e dust lanes 
were p roduced by photodissociation of H2. lAllen et al.l 
(1997) provided a qualitative analysis of similar features 
in M81 (NGC 3031), consistent with the occurrence of 
large-scale, low-density Hi PDRs connected to nearby 
bright FUV sources. They also suggested a more quan- 
titative approac h that was s u bsequ ently used on M101 
(NGC 5457) bv iSmith et all |2Q0fl), based on a simple 
equilibrium model for the physics of Hi production in 
PDRs. This provided a new method for determining den- 
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sities of molecular gas in star- forming GMCs of galaxies, 
a method that is independent of the excitation condi- 
tions for specific molecular tracers. Using a combination 
of GALEX FUV and VLA Hi data, essentially the same 
techn ique was applied to M81 (jHeiner et al.ll2008al ) and 
M83 (jHeiner et al.ll2008bt) . While evidence for the pres- 
ence of widespread molecular gas in both galaxies was 
found, we concluded that the molecular gas in the GMCs 
in M81 is generally of lower volume density. This is in 
agreement with the fact that the CO emission is gener- 
ally fainter in M81 and hence likely to be less excited 
(e.g. iKnapen et aTll2006h . 

The closer proximity of M33 allows for a much more 
detailed analysis than was possible for M81 and M83, 
and we are therefore carrying out an extensive study of 
the Hi in the immediate surroundings of OB associations 
in this galaxy in the context of the photodissociation 
model for Hi production. As an example of what can be 
accomplished with this approach, we analyse the Hi near 
two prominent OB a ssociations in M33: CPSDP Z204 
(jConrtes et alJ fl987h and NGC 604. The Hi features 
surrounding these associations are sufficiently bright in 
21-cm emission that they stand out above the general 
background, presenting a plausible case that their mor- 
phology is indeed a consequence of dissociation by FUV 
photons from the cluster of bright young stars nearby. In 
the case of CPSDP Z204 we model the Hi as a function of 
depth into a nearby GMC in terms of the "diffusion" of 
the far-UV radiation propagating into a clumpy medium. 
In all cases we determine gas densities in the surround- 
ing GMCs and compare the results to the appearance of 
CO(l-O) emission nearby where available. CPSDP Z204 
and NGC 604 have in common that local metallicity mea- 
surements are available from the literature (which im- 
proves the precision of our calculations), as well as CO 
detections (confirming the presence of molecular gas). 
We find that the higher density clouds as determined by 
our method are more likely to show CO emission nearby, 
in agreement with our earlier M83 results. 

Our method and data are described in The NGC 
604 and CPSDP Z204 findings are presented in and 
summarized in ijH 

2. METHOD AND DATA 
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We have used a simple equilibrium model of photodis- 
sociation physics to derive total hydrogen volume densi- 
ties in c andi date PDRs in M33, as we did in lHeiner et al.l 
(|2008af) and iHeiner et alj (|2008bf ). Every Hi "patch" in 
the immediate environment of the OB association is con- 
sidered a candidate PDR. The data entering into the cal- 
culation are: the measured far-UV flux of the OB associ- 
ation; the distance from the association to the Hi patch; 
the column density of the Hi patch; the local dust-to- 
gas ratio; and, some simplifying assumptions about the 
actual 3D geometry of the region. 

The following basic paramet ers of M33 were u sed: We 
assume a distance of 847 kpc (|Saha et al.ll2006f ). A UV 
foreground Galactic extinction correction of 0.33 was 
applied, appropriate for the GALEX FUV band , after 
iSchlegel et all (fl998f) and lGil de Paz et all ([20071) . 

The photodissociation-rcformation equilibrium bal- 
ance is highly sensitive to the local dust content, as it 
influences both the effective absorption cross section of 
the dust grains and the H2 formation rate. The dust-to- 
gas ratio 6 /So is normalized to the solar neighborhood 
value. We derive its value directly from the metallic- 
ity 12 + \og(0/H) a s measured in the ioni z ed gas of the 
OB assoc i ations by iRosolowskv fc Simon! ((2008), after 
llssa et aH (|1990() . It is assumed that this metallicity mea- 
surement from the H11 region is an acceptable approxi- 
mation of the metallicity in the molecular clouds that we 
are intereste d in. T his i s identical to our ap proach in 
IHeiner et all (|2008af ) and lHeiner et all (|2008bh . but now 
we can take advantage of individual metallicity measure- 
ments (where availabl e) instead of using a metallici ty gra- 
dient over the galaxy. IRosolowskv fe Simo nl lpOOl point 
out that, while their gradient fit has a high statistical 
accuracy, there is significant intrinsic scatter around this 
gradient (their Fig. 3). We will therefore assume that an 
individual measurement is a better estimate than using 
a gradient in this case. 

If we adopt a solar metallicity of 8.69 (from 
lAllende Prieto et alj 1200 If ), the dust-to-gas ratio be- 
comes: 

log 5/S = (12 + logO/H) mi - (12 + logO/H) & . (1) 

This is the same as assuming that the metallicity gradient 
is identical to the gradient in the dust-to-gas ratio, where 
the gradient is constructed from the solar metallicity (the 
dust-to-gas ratio scaled to the solar neighborhood value 
6/ So equals 1 by definition) and the value of the metal- 
licity at the location that we are interested in. However, 
since we are only interested in individual values, a radial 
dependence is not explicitel y computed. We ad opted 
NGC 604's metallicity from lEsteban et all (l2002fl and 
CPSD P Z204's metallicity from IRosolowskv fe Simonl 
(|2008h . 

The PDR model we applied comes from lAlier] (12 0041 
with improved coefficients provided by 
(2004|): 
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n(Hi) + 2n(H2) atoms cm" 3 is the total hydrogen vol- 
ume density. This gas is likely to be mostly atomic on 
the surface of t he GMC and most ly molecular deep in- 
side the cloud. lAllen et all (|2004f) (their Appendix B) 
provi de a full treatm ent of how Go is defined and de- 
rived. iHeatonl (|2009f ) recently proposed an improvement 
to this model by taking into account a more realistic 
self-shielding function of molecular hydrogen based upon 
multiple transitions instead of a single one. In that case, 
our computed volume densities need to be multiplied by 
a factor {S/5 )°' 2 , or 0.9 in the case of NGC 604 and 0.8 
in the case of CPSDP Z204, well within our estimated 
level of uncertainty. 

We used the publicly accessible GALEX FUV image 
of M33, which has a linear resolution of about 16 pc 
(4"). This data at A w 150 nm is used as a proxy for 
the dissociating UV (A « 100 nm) since the spectrum 
of the UV radiation in the ISM is fairly flat, at least 
in th e Galaxy close to the sun (jvan Dishoeck fc Black! 
The 21-cm Hi radio image of M33 was obtained 



where Vjjj is the (background subtracted) atomic hy- 
drogen column density (in cm -2 ), Go is the incident 
FUV flux measured at the affected Hi patch and n = 



from David Thilker (2007, private communication) and 
has a linear resolution of 20 pc (5"). Noise levels are 
typically 2 x 10 20 cm~ 2 . At these spatial scales we ex- 
pect to start seeing more frequently the signature PDR 
morphology in the Hi features found near FUV sources. 
The radio data were taken with the VLA and with the 
GBT (combined u sing a similar procedure as described in 
iBraun et al.1 [2009) . They are therefore photometrically 
accurate (within the noise levels). We adopt a local Hi 
background level for each region, based on measurements 
of the Hi column densities away from the candidate 
PDRs. The background levels are measured on a scale 
of typically 100 pc, and can be as high as 1 x 10 21 cm~ 2 . 
The remaining column density after subtraction of the 
background is then assumed to be produced entirely by 
photodissociating photons from the nearby OB associa- 
tion 3 . 

To apply the PDR model, we have generally followed 
the methods previously used in M81 and M83. How- 
ever, owing to the higher linear resolution, some further 
improvements have been made. The UV sources at the 
center of each region are identified (to the extent that 
they are resolved), their fluxes measured, and used to 
calculate a cumulative incident flux Go on the Hi patch 
of current interest. No interna l extinction correctio n was 
applied, for reasons outlined in lHeiner et a l. (2008a). We 
measured the distances between the central UV sources 
and a number of nearby Hi 'patches', defined as a local 
maximum in the Hi column density. These patches were 
generally selected based on their prominence in the im- 
age and/or their proximity to available CO detections. In 
these measurements we limited the distance between any 
UV source and any Hi patch to 200 pc, which is smaller 
than the cut-off of 400 pc we used previously, but more 
appropriate to the scale of the candidate PDRs we expect 
to see at this level of detail. At these separations, the 
FUV flux contributed by the young clusters of OB stars 
has also dropped to about 1% of the local ambient UV 
flux, and the additional Hi contributed by those clusters 

3 The origin of the smooth Hi background is not discussed here. 
It is possible that this too is photodissociated gas, arising for ex- 
ample from an extended population of B stars, but some fraction 
of it may be primordial. 
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Fig. 1. — Atomic hydrogen maps with selected contours. CPSDP 
Z204 (top) and NGC 604 (bottom) are shown. The central location 
of the OB associations are marked with a black cross. 

will be difficult to distinguish in the noise and confusion 
of the general Hi background. The final result of the cal- 
culation on each Hi patch is an estimate of the total hy- 
drogen volume density n = ra(Hi) + 2n(H2) atoms cm' 3 
in the GMC associated with that patch. 

3. RESULTS 
3.1. Atomic hydrogen in candidate PDRs 

As a prelude to the actual analysis of the candidate 
PDRs, we first present the atomic hydrogen images of the 
immediate surroundings of the two OB associations that 
we selected. The general morphology of the Hi will, in 
our picture, be a function of both the local surface num- 
ber density of GMCs and the ambient UV flux. Since 
we are assuming that atomic hydrogen peaks near OB 
associations are produced by photodissociation, we can 
measure these peaks and then use these results to calcu- 
late total hydrogen volume densities. It should be noted 
that we do not aim to identify all Hi patches, but rather 
attempt to find those that are clearly produced by pho- 
todissociation. 

The CPSDP Z204 area shows a variety of Hi patches 
with a suggestive shell-like feature immediately sur- 
rounding the OB association. Two more extended 
patches can be seen further away, as well as some more 
patches that either do not stand out much or are beyond 
our 200 pc cut-off. The more extended patches are stud- 
ied in more detail in this paper, because they are likely to 
be caused by the central OB association as well as being 
sufficiently resolved to try to model the Hi distribution 
radially. 



The distribution of atomic hydrogen in NGC 604 is 
part of a larger spiral arm and appears to curve around 
the central OB association in NGC 604 on the eastern 
side. While the general level of Hi column density is 
fairly high, the peaks do show a bit of structure, al- 
beit a chaotic one. We analyzed four Hi patches in this 
area, as well as considering a larger-scale approach to 
the PDR morphology by averaging the atomic hydrogen 
over a larger area. These patches were chosen partially 
because of the presence of CO in the area (a confirmation 
of the presence of molecular gas, sure to produce atomic 
hydrogen under the influence of FUV emission) . Another 
reason to look at NGC 604 in particular is the wealth of 
available data. Since our method is generally sensitive 
to low-density molecular gas, it is interesting to test it 
at the extreme end of higher densities that are known to 
occur in NGC 604 and compare the results. 

3.2. Hi morphology in the environs of CPSDP Z204 

The Hn region CPSDP Z204 (galactocentric radius 
1.63 kpc) was included in several surveys — most notably 
the o ne from which it derives its name (|Courtes et al.l 
1987), but it does not appear to have been singled out 
for individual study before. 

We have modeled the morphology of two particular Hi 
patches near FUV sources using PDR physics and sim- 
ple models of the geometry of the GMCs. We identified 
four FUV sources which define the center of Z204, and 
one additional source a bit further away to the northwest 
as shown in Fig. O We also located five Hi local max- 
ima, that we assume to be on the surface of GMCs. The 
FUV fluxes are generally about ten ti mes fainter than 
in NG C 604. CO emission detected by lEngargiola et all 
(2003) can be seen in two locations, at a resolution of 
13" or about 50 pc. The results of our measurements 
of CPSDP Z204 are listed in Table [TJ The letters and 
numbers in the finding chart correspond with the super- 
scripts of values in the table. The resulting total hydro- 
gen volume densities range from 15 to 295 cm -3 . One 
of the CO detections, close to source b, is accompanied 
by two Hi patches (nos. 4 and 5) that yield the highest 
volume densities in the area. The other CO detection, 
close to source e, does not correspond to a particular Hi 
patch at the highest column densities, but its presence 
is a confirmation that there is molecular gas in the area. 
This is important for the applicability of our method. It 
should also be noted here that the resolution of the CO 
measurements is poorer than that of the Hi map, and 
the CO location is accurate to within about 80-90 pc at 
best. 

The Hi patches in the area are consistent with being 
the photodissociated edges of GMCs, either unresolved or 
partially resolved, and illuminated by the central cluster 
of OB stars. To test this assumption we have modeled 
two patches numbered 1 and 2 in the finding chart of O 

The average Hi column densities of these patches were 
measured along the line through the center of the patch 
in the direction to the nearest cluster of UV sources. 
The column densities were averaged in thin slabs of 1" x 
12" oriented perpendicular to that line. Distances were 
deprojected in the case of patch 1 by a factor of 0.6 in 
order to obtain the separation pjjj in the plane of the 
galaxy. The resulting average column densities are shown 
in Fig. [3 
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Fig. 2. — Left panel: CPSDP Z204 at full resolution. The FUV color image is overlaid with a selection of Hi contours. These contours 
(not background-subtracted) are at 3 and 4 xlO 21 cm 2 only, and were chosen to emphasize the morphology of the Hi peaks. The noise 
level is about 2 XlO 20 cm~ 2 . Stars mark the locations where the FUV flu xes were measured. Th e orange boxes are the locations of the 
measured Hi patches. Two crosses mark the locations of CO detections by Engargi ola et all ((2003) , at a resolution of 13" (50 pc). Right 
panel: Finding chart of the same region with letters (UV sources) and numbers (Hi patches) corresponding to the results in Table [T] The 
CO detections are the crosses near sources b and e. 

somewhat "porous" to the UV radiation, which has pen- 
etrated through holes in the GMC and produced more 
atomic hydrogen further from the FUV source. A more 
elaborate version of this model will now be described for 
patch 1. 

Patch 1 looks resolved in Fig. [3J but cannot be ex- 
plained by a simple rectangular source as we have done 
for patch 2; the Hi distribution drops off too slowly 
on the right side, beyond about 100 pc (circles). The 
approximately exponential decrease at larger distances 
from the FUV source suggest a model GMC which is 
porous to FUV radiation, with many unresolved molecu- 
lar clumps embedded in a lower-density medium (see e.g. 
iFalgarone fc P hillips 199j|). These molecular clumps will 
all have photodissociated Hi on their surfaces; it is this 
Hi which we observe. The effective FUV optical depth 
of such a medium is much lower than that expected for 
uniform gas. 

We calculate first the FUV flux incident on the GMC 
at about 90 pc using the FUV luminosity of the OB as- 
sociations and the distance to the GMC. Beyond that 
point the FUV flux drops not only as the square of 
the distance, but also with an extra internal "extinc- 
tion", proportional to e~ Kr . In this case an opacity 
of K = 0.024 p c" 1 (optically thin) gives a reasonable 
fit. According to lHollenbach fc Tielen s (1999), a clumpy 
medium that is essentially empty in the interclump re- 
gions has an effective K e = pqKq. The filling factor is 
approximately po = 0.003 for GMCs that have densities 
comparable to what we are finding, which means that 
kq would be w 8; this is consistent with the properties 
of small, dense molecular clumps such as those found in 
the Orion Nebula. 

The final step in our model is to smooth this source 
distribution with the 5" beam, resulting in the solid line 
shown in Fig. [3] We note that this model provides a good 
fit to the data (circles). 

3.3. NGC 604 

NGC 604 is the largest and most luminous Hn re- 
gion in M33 and a site of active star formation, situ- 
ated at a galactocentric radius of 3. 47 kpc. It consists 
of several clusters of O and B stars (|Churchwell fc Gossl 



Fig. 3. — The average Hi column density distribution of CPSDP 
Z204 patches 1 (circles) and 2 (squares) — see the finding chart in 
Fig. [2] The dashed-dotted line represents an infinite edge fit to 
patch 1 smoothed by the telescope beam. The dotted line is the 
expected underlying distribution for photodissociated Hi produced 
on the surfaces of a clumpy GMC as the FUV radiation diffuses into 
it, and the solid line is the result of smoothing that distribution 
with the telescope beam. The dashed line is the fit to patch 2, 
assuming that this Hi patch is 3" wide in the direction of this plot 
and hence only slightly resolved by our 5" beam. See text for more 
details of the models. 



The density profile of patch 2 (squares) is the easiest 
to model. The patch seems unresolved along the line of 
sight from the UV source, and barely resolved in the per- 
pendicular direction. We fitted the data to a model of 
a simple rectangular Hi patch of constant column den- 
sity, convolved with our 5" (radio) beam. The best re- 
sult is achieved with a 3" x 9" patch with an intrinsic 
column density of 7 x 10 21 cm -2 (not including a back- 
ground level of 1.6 x 10 21 cm -2 ). The volume density 
inferred for the GMC underlying this patch is 15 cm -3 
using the observed column density enhancement of about 
3 x 10 21 cm -2 from Fig. [3l but this decreases to 4 cm -3 
after correcting the column density for beam smoothing. 
Note further that the brightness temperature associated 
with such a level of Hi column density suggests that the 
Hi is quite possibly optically thick, in which case the 
real column density may be significantly higher. The 
extra "lump" of Hi gas at about 100 pc distance from 
the FUV source might be an indication that the GMC is 
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Fig. 4. — Left panel: NGC 604 CO(l-O) image, reproduction from Tosaki ct al. (2007). The morphology of the warm dense molecular 
gas is clearly visible. Their CO(3-2) emission shows a similar morphology, tracing high excitation areas. The cross marks the center of 
NGC 604. Right panel (same scale): FUV (blue) and Hi (red) qualitative image at arbitrary scale (not background-subtracted), showing 
the distribution of atomic gas and young stars. The four boxes indicate the locations of select Hi patches. The yellow ones are close to the 
CO emission peaks. 
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Fig. 5. — NGC 604 region at full resolution. The FUV flux range is indicated in the color bar. The Hi contours are 3, 4 and 5 X 10 cm . 
The green rectangle indicates the region where the Hi columns were averaged for our measurement. The orange boxes are the locations 
of the measured Hi patches. The green crosses mark the presence of CO. See Table [l] for the measurements corresponding to the finding 
chart. 



1999), which are partially resolved in the GALEX im- 
age of NGC 604. The region is known to contain 
molecular clouds of various sizes in a clumpy m edium 
l|Lebo [ uteijlCT dTall [2001 iWilson fc Scovilld fl99l ) . Fi- 
nally, [Wilson & Scovillel (|1992l ) pointed out that M33 Hi 
data implies that NGC 604 has a relatively large mass 
fraction of atomic hydrogen, which may be the result 
of photodissociation. This makes the region an excel- 
lent target to find large-scale PDRs. We now estimate 
molecular gas densities in NGC 604 using our model, and 
compare our results with those obtained from CO maps. 

A detailed view of NGC 604 in FUV and Hi 21-cm ra- 
dio emission is shown in the right panel of Fig. [U The 
FUV emission (blue) is spread across the region, but it 
is strongest at the Hn region defining NGC 604. The 
atomic hydrogen (red) follows one of M33's spiral arms, 
and appears to "wrap around" to the east of NGC 604, 
forming a large Hi structure of typical size 200 pc similar 



to tho se structures first identified in M81 by I Allen et al.l 
(|1997f l. The Hi to the west of the arm drops below the 
detection limit, but it is present at faint levels in lower 
resolution maps. Five resolved OB associations were 
identified and their FUV fluxes measured. Four repre- 
sentative Hi patches were selected as candidate PDRs, 
two of them based on the detection of CO close to their 
locations. This is shown in Fig. [5] Table[T]lists the mea- 
surements of candidate PDRs, taken in a similar fashion 
as in CPSDP Z204. The resulting total hydrogen den- 
sities of the underlying GMCs vary from about 20 to 
140 cm~ 3 , with the highest values corr esponding to the 
two lo cations near the CO detections bvlEngargiol a et ail 
(|20Q3l) . 

We have also analyzed the Hi "envelope" around NGC 
604 as if we had observed it's Hi with a typical resolution 
of ~ 100 — 200 pc, analogous to our earlier work on M101, 
M81 and M83 (boldface in Table Cp). At this linear res- 
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olution, the (smoothed) Hi column densities are lower, 
leading to higher GMC volume densities of typically 500 
cm~ 3 estimated with a simple application of the PDR 
model. However, the more detailed analyses we have 
provided in this paper show that the true mean densities 
are likely to be significantly lower when a clumpy source 
geometry and the effects of finite telescope resolution are 
included 4 . 

The densities of the GMCs that we found here are lower 
than those obtained using CO methods. This is caused 
by our method being specifically sensitive to lower densi- 
ties. Higher density clouds have a thinner layer of atomic 
hydrogen, as well as a lower beam filling factor as they are 
more concentrated. The GMCs with the highest densi- 
ties that can be detected through their CO emission will 
therefore correspond to smoothed Hi column densities 
that are too low to stand out from the general back- 
ground level. What can be seen here is that our method 
yields cloud densities on the high end of the range of 
densities to which we are sensitive. 

4. CONCLUSIONS 

1. Using VLA + GBT observations of the 21-cm Hi 
line and GALEX observations of the FUV emission, 
along with a simple model for the photodissociation 
of molecular hydrogen on the surfaces of GMCs, 
we have estimated total hydrogen volume densities 
in the GMCs surrounding the two OB associations 
NGC 604 and CPSDP Z204 in M33. 

2. We find values for the molecular gas densities in the 
parent GMCs near these OB associations ranging 
from 15 to 540 cm~ 3 . Regions with higher gas den- 

4 This apparently counter-intuitive behavior is a consequence 
of Equation [2] which shows that for a given Go the only way to 



sities (> 100 cm 3 ) appear to be correlated with 
the presence of CO(1-0) emission. 

3. We explain the morphology of two Hi patches in 
the environs of CPSDP Z204 using straightforward 
assumptions about the underlying geometry of the 
molecular gas. One patch appears to be a slightly- 
resolved GMC. The other source was largely re- 
solved; it can be modeled as a clumpy molecular 
medium which is very porous, and therefore sig- 
nificantly more transparent to FUV emission than 
would otherwise be expected. In our model, the Hi 
we are observing is produced by photodissociation 
of H2 on the surfaces of these molecular clumps. 

Our approach views the Hi found in the environs of 
star-forming regions of galaxies to be a product of the 
star formation process, not a precursor to it. The ability 
to model the morphology of Hi features in the environs 
of OB associations and to predict the appearance of CO 
emission as an excitation effect related to higher volume 
densities in the parent GMCs are two results of this pa- 
per which lend further support to this approach. A more 
extensive treatment of GMCs in M33 based on this ap- 
proach is in preparation (Heiner 2009, PhD thesis). 
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increase iVjjj is to decrease the total volume density n. 
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TABLE 1 

Measurements of NGC 604 and CPSDP Z204 



NGC 604 CPSDP Z204 



<5/<5o 0.63 0.38 

FUV fluxes (10~ 15 ergs cm" 2 s" 1 A" 1 ) 406°, 15. 3 b , 22. 6 C , 3.47 d , 914 6.03°, 6.73 b , 9.61 c , 3.72 d , 10.9 e 

8.54 e 

N bg (10 21 cm" 2 ) 1.0 1.0 

"^H (10^^ Cm - ^) 4 2^ :a ' ,cc ^ 2 f}2:ade ^ Q3:abcd ^ g4:afecd 2 4 Q^ :a ^ C( ^ 4 4.2 :e 3 g3:abcd 

(background-subtracted) 3 A i:abcde , 3.2 5:a6cd 

Go (cumulative) 7.40 1 , 6.06 2 , 10.53 3 , 6.93 4 , 15.21 0.87 1 , 0.65 2 , 0.45 3 , 7.20 4 , 3.97 5 
G/G bg range 0.01 - 4.45 1 ; < 0.01 - 4.01 2 ; 0.01 - 6.27 3 ; 0.02 - 0.18 1 , 0.57 2 , 0.01 - 0.09 3 , 

< 0.01 - 4.45 4 ; 6.49 0.13 - 1.65 4 , 0.05 - 0.51 5 

n (derived, in cm" 3 ) 34 1 , 109 2 , 139 3 , 22 4 , 539 18 1 , 15 2 , 16 3 , 295 4 , 180 5 

Fractional error range 0.26 - 0.44 0.26 - 0.36 



A description of this table can be found in i) £|3.2l and l3.3l 



